HUTP-06/A0037 
LBNL-61548 
MCTP-06-20 



Disentangling Dimension Six Operators 
through Di-Higgs Boson Production 

Aaron Pierce 1,2 , Jesse Thaler 3 ' 4 , Lian-Tao Wang 1 

1 Jefferson Laboratory, Physics Department, Harvard University, 

Cambridge, MA 02138 

2 Randall Laboratory, Physics Department, University of Michigan, 

Ann Arbor, MI 48109 

Department of Physics, University of California, Berkeley, CA 94720 

4 Theoretical Physics Group, Lawrence Berkeley National Laboratory, 

Berkeley, CA 94720 



Abstract 

New physics near the TeV scale can generate dimension-six op- 
erators that modify the production rate and branching ratios of the 
Higgs boson. Here, we show how Higgs boson pair production can yield 
complementary information on dimension-six operators involving the 
gluon field strength. For example, the invariant mass distribution of 
the Higgs boson pair can show the extent to which the masses of ex- 
otic TeV-scale quarks come from electroweak symmetry breaking. We 
discuss both the current Tevatron bounds on these operators and the 
most promising LHC measurement channels for two different Higgs 
masses: 120 GeV and 180 GeV. 



1 Introduction 



If the Large Hadron Collider (LHC) does not produce any resonances aside 
from the Higgs boson required to unitarize W-W scattering, physicists will 
be forced to look for new physics in indirect ways. One approach, recently 
re-emphasized by [T] , is to hunt for new physics via the presence of higher- 
dimension operators involving only Standard Model fields. Many of these 
operators, exhaustively catalogued in [2], are already well constrained by 
existing precision measurements from LEP and are unlikely to be probed 
further at the LHC. Here we discuss higher-dimension operators containing 
the Higgs boson that are currently poorly constrained, but could directly 
influence collider phenomenology at the LHC. Our primary focus will be on 
final states with two Higgs bosons. 

Colored particles that get part of their mass from electroweak symmetry 
breaking (EWSB) can induce the operator 

O^c^G^rfH (1) 

at the loop level. The mass scale v = 246 GeV has been chosen for later 
convenience, and Attv may or may not be the actual scale of new physics. The 
influence of this and other operators on single Higgs boson production and 
branching ratios was recently discussed in [H [3] . By itself, 0\ is insufficient 
to completely describe the low energy effects on both single and pair Higgs 
boson production. To see this, consider a new particle whose mass comes 
entirely from EWSB. This yields a different (non-decoupling) operator. As 
is familiar from Higgs low energy theorems, a heavy quark with Yukawa 
coupling A — > oo generates not 0\ but 



a 



(H*H\ 



°* = Pa l°g (^r J, (2) 



which can be understood by thinking of if as a background field and treating 
the heavy quark mass as a threshold for the running of the QCD gauge 
coupling |4|. If we expand 0\ and O2 in terms of the physical Higgs boson 
h(H = ±(h + v)), 

0i ,S£g^Q + *), 0,^0^(1-*), (3) 
then the differing effects on Higgs boson pair production are manifest. 
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It is also clear that these two operators are sufficient to describe single 
and pair Higgs production at energies well below the mass scale of the new 
physics. General models are effectively described by a linear combination 
of 0\ and Oi- A study of Higgs boson pair production would indicate the 
relative importance of these two operators. Such an observation would probe 
the extent to which new colored particles receive their mass from EWSB. 

There is even a parametric limit in which deviations arising from the oper- 
ator 0\ might become visible before any new particles are directly observed. 
Imagine that there exist Nf new colored particles, all with heavy masses and 
some coupling to the Higgs boson. Then the cross section for their direct 
production scales like Nf, but processes that involve the operator 0\ will 
go like Nf. So, in the limit of large Nf, the indirect effects of the operator 
will be visible first, and the effects described in this letter could be the first 
indications of new physics. 

Of course, even if the new physics is produced directly, Higgs pair produc- 
tion would remain an interesting channel to help disentangle the new physics. 
New physics predictions for Higgs pair production already exist in a model- 
dependent context, including Little Higgs Models [5], Randall- Sundrum like 
models [6j, extended Higgs Sectors [HE], top condensation models [9], and 
the Minimal Supersymmetric Standard Model [10J [EE]- By using the lan- 
guage of 0\ and O2 we can show the importance of Higgs pair production 
mo del-indep endent ly. 

In the following section, we describe the new physics that can give rise 
to the operators 0\ and Oi- Next, we review the constraints on the possible 
size of these operators. We then discuss possible discovery of these operators 
in Higgs pair production, considering two concrete cases: rrih ~ 120 GeV 
and m/i ~ 180 GeV. We conclude with a brief discussion of other higher- 
dimensional operators involving the Higgs boson. We find that the operators 
considered here are the ones most likely to give rise to interesting effects at 
the LHC. 

2 Theoretical Considerations 

The operators 0\ and O2 have a different dependence on m? hh than the Stan- 
dard Model contribution to Higgs pair production. In particular, amplitudes 
involving will grow like m\ h all the way up to the mass scale of the 
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Figure 1: The contributions to Standard Model Higgs pair production are 
dominated by loops containing top quarks. 

new physicsfj] On the other hand, contributions to the amplitude from Stan- 
dard Model processes shown in Fig. [I] do not grow for energies above the 
top quark mass. Assuming the scale of the new physics is much larger than 
m t , this dependence on m 2 hh will serve as a way to disentangle the effects of 
higher- dimension operators from the Standard Model contribution. 

As shown in Eq. ([3]), 0\ and O2 induce different contributions to single 
and pair Higgs boson production. Thus, different kinds of ultraviolet physics 
can yield different relative signs between the operators containing one and 
two physical Higgs bosons. In general, there will be interference between the 
two diagrams in Fig. [3 The amount of interference will give us a handle 
on the relative size of c\ and c-i. A similar interference is well-known in the 
Standard Model [12] and could potentially be used to measure the h 3 Higgs 
self-coupling [13], [H] . 

There is no renormalizable Lagrangian that can generate 0\ at tree level. 
To see how 0\ and 2 can be generated at loop level, consider the effect of 
new heavy quarks who get some of their mass from EWSB. If we treat H 
as a background field, then each quark mass is a threshold for the one-loop 
running of the QCD gauge coupling [I]. Assuming all of the heavy quark 
masses rrii 3> m^, the low-energy QCD gauge coupling is 

1 1 buy A 1 , rmdH)) 



g 2 (fi) g 2 (A) 8vr 2 & /x 8vr 2 ^ 1 & /1 



1 In this paper, we will assume that the operators 0\ and O2 completely dominate the 
new physics contribution at the LHC, and that higher derivative operators that might 
soften this behavior are subdominant. 
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Figure 2: The two diagrams that contribute to Higgs pair production coming 
from the higher dimension operators 0\ and 02 • In the first diagram, a new 
g-g-h vertex combines with the Standard Model three Higgs boson coupling. 



where Sbj = 2/3 for a SU(3) fundamental fermion. The non-canonically 
normalized QCD gauge kinetic term 

^kinetic = 4g2^^M^a (^) 

can then be expanded in terms of h to determine the effective values of C\ 
and C2 relevant for single and pair Higgs production!! At higher order in h, 
0\ and O2 are insufficient to specify all of the allowed Higgs interactions, so 
Eqs. fll]) and (jSJ) should be used directly. 

For concreteness, consider the following Lagrangian (H = e • W) 

-£mass = Ai (QHT C + QHB C \ + A 2 [Q C HT + Q C HB) 

+ m A QQ c + m B (TT c + BB c )+h.c, (6) 

where Q, Q c are vector-like SU{2) L doublets and T, T c (B, B c ) are vector-like 
SU{2) L singlets, with appropriate hypercharges and SU(3)c couplings. In 
order to suppress contributions to the T-parameter [T5] , we assume custodial 



2 Note that the definitions of 0\ and O2 assume canonically normalized kinetic terms. 
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isospin. Using Eq. (jl]) with i = 1 to 4 and <56j = 2/3 El 



4-/5 4 1 



= 



2mATriB 



(7) 



3 (l-/?) 2 ' ° 2 3 (l-/?) 2 ' 



AiA 2 f 2 ' 



If all the mass of the heavy quarks comes from EWSB ((3 = 0) then c\ = 0. 

Can the effects of 0\ be visible before the new colored states are seen 
directly? In the case of heavy quarks that get all of their mass from EWSB, 
it seems unlikely. The new quarks could at most have Yukawa coupling 
A ~ 47r / a/ Nc to keep the theory perturbative, where the number of colors 
N c = 3 for our toy model. With masses of \v/y2 ~ 1.3 TeV, the heavy 
quarks will have a rather small pair-production cross section ~ 10 fb, but 
could well be visible at the LHC in single production via b — W fusion (see, 
e.g. [16J), depending on the flavor structure of the heavy sector. So, direct 
production would likely be the first window on new physics of this type. What 
about the case where the quarks have mostly vector-like masses {(3 ^> 1), 
so Ci ^> C2? The large vector-like mass supresses the overall contribution 
to ci, as this operator decouples like m 2 . So for c\ to be 0(1), the large 
vector-like mass must be compensated by a large number Nf of heavy quarks: 
Nf ~ m 2 /X 2 v 2 to prevent rapid decoupling. Since the b — W fusion process 
scales roughly as m -7 , the total production cross section for Nf copies of new 
physics will scale like m -5 . Thus, there is at least a parametric limit where 
the effect of 0\ is visible before the new heavy quarks are seen directly. For 
reasonable values of A, (3, and Nf, the mass of the new quarks could even 
exceed the LHC center-of-mass energy while the contribution of these heavy 
states to 0\ could remain substantial. 



Direct experimental constraints on 0\ and O2 are quite weak for a low mass 
Higgs. Direct constraints come from Tevatron searches for Higgs boson pro- 
duction via gluon-gluon fusion, which constrain the combination (ci + C2). 

3 In order for this QCD beta function argument to make sense, f3 has to be far from 1, 
or else there is a mass eigenstate lighter than the physical Higgs. Note that (1 — /3) 2 is 
proportional to the determinant of the mass matrix. For simplicity, we take (3 as a real 
parameter. In the case where it is complex, the formulae for c\ and ci are modified, but 
they remain real, as Hermiticity of the Lagrangian requires. When there are phases in the 
mass matrix, ci can be negative. 



3 Experimental Constraints 
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This production channel is generally not useful at the Tevatron when search- 
ing for Standard Model Higgs boson with unmodified properties — the back- 
grounds are generically too large. In the low- mass (rrih ^ 130 GeV) re- 
gion, the Higgs boson dominantly decays to b quarks, with the leading sub- 
dominant decay to r leptons. Decays to b quarks cannot be used for the 
Higgs boson search because of the too-large QCD background from di-jets. 

A more promising channel is h — > 77. In the Standard Model, this 
branching ratio is small, Br(/i — > 77) ~ 2 x 10~ 3 . Current bounds from 
the DZero experiment constrain the branching ratio to be less than Br(/i — > 
77) < 0.5 [17], assuming the Standard Model production cross section. To 
see how this measurement constrains (07 + C2), it is convenient to write the 
modified cross section times branching ratio for single Higgs production as 

a(i - h) x Br(fc - f) = aSM{ ^ k) ^l, (8) 

i 

where Tij is the partial width for Higgs decay into the i and / states, and F 
is the total width: 

r + r 77 ~\~ T ww + T Z z + Tbfc + r cc + Trr + ■ ■ • • (9) 
In this language, the DZero constraint can be written as 

T " < 0.5. (10) 



■PSM 
99 

If the partial width T 77 is unchanged from the Standard Model, then Eq. (fTUl) 
implies no bound on T gg , and therefore no bound on (ci + c 2 ). 

In actuality, we expect the new physics that contributes to T gg to also 
modify T 77 through operators like 

11 = c,^F^h. (11) 

Because T gg starts to dominate the width when (C1+C2) < —1.75 or (C1+C2) > 
1.05, for large enough \c\ — Cz\, the branching fraction to photons will be 
controlled by the ratio of the color charge to the electric charge of these new 
statesQ For particles with top and bottom quark-like quantum numbers, we 



4 The Standard Model contribution to T 77 from the Vl^-loop will be negligible in this 



regime. 
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expect the asymptotic behavior 
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77 
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(12) 
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and Eq. ffTUl) gives the bound: 

— 2.8^, (ci + c 2 ) ^, 2.1, (Direct Search for gg — > /i — ► 77). (13) 

To obtain these values, we have taken into account the momentum depen- 
dence of the top quark loop (a 6% effect for a 120 GeV Higgs), and assumed 
similar NLO K-factors for the Standard Model contribution and the contri- 
bution coming from the new operator. This is likely a good approximation, 
as the radiative corrections are dominated by interactions involving the ini- 
tial state gluons. The asymmetric bounds on (ci + c 2 ) show the effect of 
constructive versus destructive interference with the Standard Model. 

In passing, we note that the search for Higgs bosons in the h — > tt 
channel is unlikely to add any useful new constraints on this operator. For 
a rrih ~ 120 GeV, searches for the A of the MSSM to a pair of r leptons 
places a limit a(pp — ► h) x Br(/i — ► t + t~) < 15 pb [18J. The rate for this 
process in the Standard Model is about 0.056 pb. The decay rate h — ► t + t~ 
is dominated by SM tree level process and not sensitive to the class of new 
physics we consider here. Therefore, we do not expect any constraint on 
(ci + c 2 ) from Tevatron search of this process. A similar story applies to 
searches for gg — + h — > WW* in this mass region: by the time the production 
cross section is large enough to be observable, the Higgs width is dominated 
by gluons, and the branching fraction to W's is too small to be observed. 
Also, even accounting for possible branching ratio enhancements, the search 
channel h — > 7Z gives slightly worse bounds than the diphoton channel (19] . 

When the Higgs is heavier than 160 GeV, the bounds on (ci + c 2 ) be- 
come more severe [20]. Even accounting for the possible increase of T gg , the 
dominant decay mode is to a pair of on-shell W's at tree level. Since the 
total width is dominated by IVvk? from Eq. (jSj) we see that the rate of this 
process is proportional to (ci + c 2 ). The combined bounds from the Tevatron 
allow an increase in the Higgs production cross section of a factor of 6 over 
the Standard Model rate, corresponding to 



- 1.2 < (ci + c 2 )< 0.5, (Search for gg -> h -> WW, m h ~ 180 GeV). 



(14) 
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Similar bounds apply throughout region > 160 GeV. For this size of Ci, 
the partial width Y gg does not compete with Ywwi so the Higgs branching 
ratio to W bosons is largely unchanged. 

While the above discussion summarizes the state of the direct bounds on 
(ci + C2), there are potential limits from precision electroweak measurements. 
In principle, the presence of 0\ and O2 need not generate operators at a 
dangerous level, but it should be noted that heavy states that contribute to 
0\ and O2 could contribute to S and T. Even if we assume custodial isospin, 
new particles can contribute to the operator 

Owe = ^W^HB^, (15) 

which is directly related to the S'-parameter [15]. The exact contribution to 
Owb will depend on the electroweak quantum numbers of the new particles. 
If the new particles have quantum numbers identical to Standard Model 
quarks, this bound on S will give a fairly strict limit on the potential size 
of 0\ or Currently, electroweak precision measurements constrain S = 
—0.13 ± 0.10 [21]. The contribution of a fourth generation of quarks to Owb 
yields 

AS=±. (16) 

Rescaling the quark masses to force AS" < .10 implies the following rough 
bound: 

AS < .10 Ci < 0.4. (17) 

Notice that there is no parametric limit in which we could both have sizable 
contributions to 0\ t 2 and decouple contributions to Owb- For example, 
in our toy model, Nf will scale like m 2 /\ 2 v 2 , which cancels the high scale 
suppression in front of the H^W^ V HB^ U operator. 

However, if the new colored states have exotic quantum numbers, there 
can be a sizable enhancement to the c$. In particular, the contributions to q 
goes like the Dynkin index of the SU(3) representation, while the contribu- 
tion to the S-parameter goes like the the dimension of the SU(3) representa- 
tion. For exotic representations, this can cause a substantial deviation from 
the prediction of Eq. (1171) El Furthermore, it is possible that exotic representa- 

5 Sufficiently large representations will cause SU(3)c to become asymptotically non- 
free, and hit a Landau pole slightly above the mass of the exotic representation. We 
remain agnostic as to what new physics lies above the scale where these new operators are 
generated, so we do not view this constraint as too limiting. 
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tions of SU(2) L x SU (2) R might be present, and these can give contributions 
to S of either sign [22J, so a combination of representations might leave S 
unchanged while giving a large contribution to 0\ and 2 . Finally, colored 
scalars can contribute to 0\ without contributing at all to the S parameter: 
the interaction ^(pWH does not require cf) to have any electroweak quantum 
numbers. So, while q ~ 0.4 may represent a typical value in some theories, 
values much larger are certainly possible, so we will consider the q as free 
parameters up to the limits imposed by the direct searches, Eqs. (TL3]) and 



4 Measuring Higgs Pair Production 

The sizes of c\ and c 2 could be determined uniquely by measuring the cross 
sections cr(pp — > h) and cr(pp — > hh). However, the quantity that is most 
directly measured experimentally is not the cross-section, but rather a cross- 
section times a branching ratio. Given the likely possibility that new physics 
will modify the Higgs branching ratios, it is useful to have an independent 
measure of c\ and c 2 . 

A differential distribution of the form 



where x is some kinematical variable, such as m\ h , gives such a handle. In 
general, f(c\, c 2 ) is a function that depends on the size of ci, c 2 , their relative 
sign, as well as interference with the Standard Model piece. Therefore, both 
the rate and the shape of the di-Higgs distribution give independent probes 
of the coefficients C\ and c 2 , with different systematic errors associated with 
each measurement. 

Whether the shape of the di-Higgs invariant mass distribution can be 
observed or not depends on the overall rate of Higgs pair production at a 
given luminosity as well as the Standard Model background to the channel 
used to reconstruct the Higgses. In the following subsections, we comment 
on two different Higgs mass windows: a low-mass Higgs boson near 120 GeV 
and a Higgs boson above the W + W~ threshold near 180 GeV. To investigate 
this question, we augmented MadGraph [23] with new HELAS [21] routines 
to simulate the contributions of both the Standard Model top quark loop 



da 

dx 



f(0l, 02) X ^ X ^, 



(18) 
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H21 123 E3 and the operators O x and 2 E The FF package [2B] was used 
in the numerical evaluation of the relevant top loop integrals. In the limit 
where the contributions of our new operators vanish, our numerical results 
agree with those of [10J. 

A search for 0\ and O2 would involve the same final states as the search 
for the h 3 Higgs self coupling [T31 EE], but typical values of C\ and c 2 can 
lead to an order of magnitude increase in the cross section cr(gg —> hh) (and 
more extreme values of q can give enhancements of even larger factors that 
are nevertheless consistent with all known phenomenological constraints). A 
preliminary ATLAS [27] study suggests that the Higgs self coupling could 
only be measured with luminosities typically associated with the SLHC It 
would be interesting to know whether the possible drastic increase in cross 
section from dimension six operators would paint a more optimistic picture 
for Higgs pair production at the LHC even after accounting for detector 
resolution effects. 

4.1 Low Mass Region 

In the region where the Higgs has a low mass ~ 120 GeV, two Higgs 
decays can best be observed via the process gg — > hh — > 6677 where at 
least one of the jets is 6-tagged [28]. While this process has manageable 
background, in the Standard Model it suffers from the small branching ratio 
Br(/i — > 77) ~ 2 x 10~ 3 . As discussed in the previous section, both the 
branching ratio to photons and overall production rate can be substantially 
changed from the Standard Model prediction. However, the branching ratio 
to 6-quarks will generically decrease because of the increased partial width 
T gg , and at the extremes of the Tevatron allowed region in Eq. (fT5|) . Br(A — > 
bb) decreases by about factor of four. Together with a possible factor of five 
increase in Br(/j — > 77) described in Eq. (fl2l . the overall branching ratio of 
hh — > 6677 could be comparable to the Standard Model value. Away from the 
extreme values of q, it is possible that the product Br(/i — > gg) x Br(/i — > 77) 
could be larger or smaller than the Standard Model, depending on the extent 
of cancellation between the new physics and the iy-loop contribution to the 
h — > 77 rate. 

The enhancement of a(gg — > hh) relative to the Standard Model value for 

6 We thank Rikkert Frederix for providing the HELAS routines to implement the 
G^G^h vertex. 
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Figure 3: The ratio of cr(gg — > hh) to the Standard Model di-Higgs cross 
section for rrih = 120 GeV. This includes the effect of interference between 
the contributions from Oi, O2, and the Standard Model. We assume the new 
contribution to di-Higgs production inherit the same NLO K-factors as the 
Standard Model. The Standard Model cross section is 30 fb, and the allowed 
range in Eq. (TTBl) from direct Tevatron searches is —2.8^ (c\ + c 2 )^ 2.1. 

a range of values c\ and c 2 is given in Fig. [31 Each tree-level cross section is 
multiplied by a K-factor of 1.65 to take into account NLO effects [21], where 
we are assuming that the QCD corrections to the diagrams in Figs. [T] and 
[2] are comparable. The factorization scale and the renormalization scale are 
taken to be rrihh, and the CTEQ6L1 parton distribution functions [30J are 
used. As has been previously noted, the dependence on the renormalization 
scale is somewhat severe for this case, and can easily affect the cross section 
by tens of percent. 

In [28], it was found that with a set of reasonable cuts to isolate the signal, 
there could be 6 signal events compared to 14 background events with 600 
fb -1 of LHC data. Given the large possible enhancements of the rate of 
Higgs pair production in Fig. [31 the signal to background ratio can be much 
enhanced, and the event rate need not be so tiny. This assumes that 0\ and 
O2 do not make large contributions to the background. In particular, one 
might worry that an enhanced g—g—h vertex could influence the irreducible 
background from hbb and /177. However, as shown in [28J, these backgrounds 
are negligible compared to the background from QCD with fakes, so it is 
reasonable to expect that the values in Fig. [3] represent the real increase in 
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the signal to background ratio except at the most extreme values of the c,. 

The differential cross sections as a function of are shown for various 
choices of the Cj in Figs. H] and [51 Because the rate of single Higgs production 
is related to (C1+C2), we keep the sum (01+02) fixed in Fig. [5]while varying the 
difference. The large variety of shapes and normalizations show how di-Higgs 
production can be an important probe of the properties of the Higgs boson. 
There are two specific features to notice in these figures. First, compared to 
the Standard Model alone, there is a long tail in the rrihh distribution when 
Oi is turned on, and the size of that tail is governed by \c\ — C2 1 - Second, 
for values of C\ and C2 that are not too large (< 0.5), interference with the 
Standard Model can be important, possibly causing a significant deficit of 
events around m hh = 400 GeV. 

With sufficient luminosity, a measurement of cr(gg — ► hh) as well as of 
the mhh distribution in di-Higgs events will give some handle on the sizes 
of 0\ and Oi- It is beyond the scope of this paper to try to estimate the 
errors on the measurements of c\ and C2 due to energy resolution, background 
subtraction, and statistics. However, because of the large enhancement of 
di-Higgs production due to these new operators, it is not unreasonable to 
expect some m^h shape information to be available after several years of 
high luminosity (100 fb _1 /yr) running at the LHC. 

4.2 Above the W+W Threshold 

For a Higgs with a mass near 180 GeV, its dominant decay mode is to W + W~ . 
The hh —> W + W~W + W~ mode was considered in [T3l [T4] . where it was 
determined that the the cleanest channel for discovering Higgs pairs was when 
the four W*s yielded two same sign leptons, i.e. hh —>■ W + W~W + W~ — > 
jjjji v^u. For various values of Cj, the cross sections in this channel are 
given in Fig. [6j Using the cuts in [H] with 600 fb _1 , there are 110 Standard 
Model di-Higgs events with an expected signal/background ratio of 1 to 5. 

Because of the presence of two final-state neutrinos, the invariant mass 
of the two Higgses cannot be determined uniquely. In [H], an approximate 
variable m visi bi e was used, which systematically underestimates the real in- 
variant mass by taking the invariant mass of just the final state leptons and 
jets. On the other hand, once jets are combined to give pseudo-H^'s, the 
decay topology of hh — > W + W + £~b>£~v allows for a full reconstruction of 
the neutrino four-vectors, up to discrete ambiguities. Eight constraints are 
required to measure the eight components of the two neutrino four-vectors. 
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Figure 4: Differential cross-sections as a function of m^h- In the top graph 
C2 = is fixed while c\ varies, and in the bottom graph c\ = is fixed while 
c 2 varies. We have set m h = 120 GeV and m t = 174.3 GeV. The curves 
for rrih = 180 GeV are quite similar, with the trivial modification that the 
threshold energy is changed. Note that the asymptotic behavior at large 
is controlled by the difference \c\ — c%\. When c\ = —0.5 and c 2 = 0, there is 
a pronounced dip at = 400 GeV, coming from interference between 0\ 
and the Standard Model top loops. 
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Figure 5: Differential cross-sections as a function of rrihh for nih = 120 GeV. 
Here, the single Higgs production rate is fixed by fixing (c\ + c 2 ), but the 
properties of di-Higgs production are clearly modified as the proportion of 
0\ and O2 changes. 
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Figure 6: The ratio of cr(gg — > hh) to the Standard Model di-Higgs cross 
section for nth = 180 GeV. The Standard Model cross section is 13 fb, and 
the bounds in Eq. ([14"]) from direct Tevatron searches is — 1.2< (C1 + C2) ^ 0.5. 
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They are obtained by requiring that the two neutrinos are on-shell (2), that 
the neutrinos and charged leptons must reconstruct two W's (2), that oppo- 
site sign W pairs must reconstruct two Higgses (2), and that the transverse 
momenta of the neutrinos must yield the missing px vector (2). 

Of course, finite width effects, initial and final state radiation, and energy 
resolution issues will blur the di-Higgs invariant mass distribution. Even with 
a perfect understanding of those experimental issues, there can be anywhere 
from zero to eight real solutions to the constraint equations^ and one must 
pick a method to choose the right neutrino four-vectors. Here, we consider an 
ad hoc method to determine rrihh that seems to give reasonable results with 
high efficiency: each solution to the constraint equations yields a different 
value of rrihh, and we plot the mean value of rrihh in events where the standard 
deviation of rrihh is less than 10% of the mean. This technique removes around 
half of the candidate events, and ensures that only values close to the true 
values are plotted. As shown in Fig. this reconstruction method does not 
appear to systematically distort the rrihh distribution and gives a much better 
estimate of invariant mass distribution compared to the m V i S ibi e variable. 

One of the important benefits of this reconstruction technique could be 
on background reduction. One of the largest sources of background comes 
from pp — > Whjj where the two jets appear to reconstruct a W. In the signal 
sample, the failure rate for the neutrino reconstruction (i.e. when there are 
zero solutions to the constraint equations) is around 20%. However, in a 
background sample of pp — > Whjj where the invariant mass of the two jets 
(ignoring smearing) are forced to lie within 10 GeV of the W mass, the 
failure rate is nearly 80%. This is likely because the largest contribution to 
pp —>■ Whjj comes from diagrams involving Higgs-strahlung off of W's, which 
has a very different kinematic structure from true Higgs pair production. It 
would be interesting to see whether this four-fold improvement of the signal 
to background ratio persists in a more realistic background study. 

As mentioned in Figs. H] and [5j the theoretical rrihh distribution at rrih = 
180 GeV is qualitatively similar to the rrih = 120 GeV case apart from 
the different kinematic threshold. It is beyond the scope of this paper to 
consider the full effect of energy resolution and QCD radiation on Higgs pair 

7 Zero solutions arise if enough of the Higgses and W's are sufficiently off-shell. In this 
case, the kinematics of the event will not be consistent with, say, the central value of the W 
mass of 80.3 GeV. To get eight solutions, note that both the W and Higgs reconstruction 
equations are quadratic, yielding a maximum of four real solutions, and there is a two-fold 
ambiguity as to which hadronic W should be paired with which leptonic W. 
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Figure 7: Distributions of rrihh at rrih = 180 GeV for (top) c\ = 0, C2 = 0.5 
and (bottom) c\ = —0.5, C2 = 0. The neutrino reconstruction technique does 
a better job at matching the real invariant mass distribution than the risible 
variable, at the expense of reducing the statistics by a factor of two because 
of the 20% failure rate for finding any solution to the constraint equations 
as well as our ad hoc method for resolving reconstruction ambiguities. In 
these plots, the curves are normalized to peak at 1 in arbitrary units. Note 
that regardless of whether one uses the visible or reconstructed mass, the 
shape difference between the top and bottom curves is still observable, so 
the main reason to use the reconstruction technique is to potentially reduce 
backgrounds. 
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reconstruction. Because of the large number of final states and the need to 
accurately know the missing px vector, neutrino reconstruction may or may 
not be the best experimental technique for extracting information about the 
di-Higgs rrihh distribution. Most likely, the most appropriate technique will 
depend on the total number of events seen. 

5 Conclusions 

If new colored resonances are observed at the LHC, the signature explored in 
this paper will be useful to explore the extent to which these new particles get 
their mass from EWSB. If we are unlucky at the LHC and new particles are 
inaccessible directly but just around the corner, then the di-Higgs invariant 
mass distribution can be used to probe them indirectly, possibly providing 
the first measurement of physics beyond the Standard Model. 

One might wonder whether there are other operators that could con- 
tribute to Higgs pair production at the LHC. For example, consider the 
following dimension six operators from [2J: 

O fn = iPH X c H^H, O z > = ipa^H^D^H, (19) 

where if) represents a Standard Model SU{2) L doublet, and \ c a singlet. We 
might expect to generate Ofn in some Froggatt-Nielsen [HT] model, and Oz> 
could arise from integrating out a Z' gauge boson that coupled to Stan- 
dard Model fermions and the Higgs. However, both of these operators are 
constrained by existing new physics searches. Unless we assume minimum 
flavor violation, then Ofn can introduce dangerous flavor-changing neutral 
currents, but minimal flavor violation implies that the coefficient of this op- 
erator will be tiny for the first two generations. Similarly, setting the Higgs 
to its vacuum expectation value in Oz> will generate anomalous couplings 
between Standard Model fermions and the Z, which are well-constrained by 
LEP. 

At best, we could consider Ofn applying only to the third generation. If 
we tune a modified renormalizable Yukawa coupling qHb c against the opera- 
tor qHb c H^H, we can keep the bottom mass fixed while generating large bbh 
and bbhh couplings. While one could study the interference between these 
two couplings in Higgs pair production, we know of no well-motivated model 
to justify the necessary tuning. For a more motivated scenario, the Yukawa 
coupling qHb c could be eliminated altogether, and the bottom mass could 
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come solely from Ofn- This is natural in a two Higgs doublet model, and 
the smallness of the bottom mass relative to the top would be explained by 
the fact that symmetries force the bottom coupling to the Higgs to come 
from a dimension six operator. However, in this scenario with tan (3 = 1, the 
bb — > hh cross section at the LHC is on the order of 10 ab, too small to be 
relevant. 

Thus, Oi and 2 are selected out as particularly interesting operators: 
these new contributions to TeV-scale physics are not currently constrained 
by experimental searches, but they have the potential to induce novel physics 
in Higgs boson pair production at the LHC. Even if the LHC experiments 
do not find any new resonances beyond the Higgs boson, high luminosity 
studies of the Higgs properties at the LHC could still offer a glimpse of the 
ultraviolet. 
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